Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 4 May 2009 (MN MfcjX style file v2.2) 



Starburst— driven galactic outflows 



o 
o 



O 

u 
6 



Biman B. Nath 1 & Joseph Silk 2 

1. Raman Research Institute, Sadashiva Nagar, Bangalore 560080, India 

2. Department of Astrophysics, University of Oxford and Beecroft Institute for Particle Astrophysics and Cosmology 
Denys Wilkinson Building, Keble Road, 0X1 3RH Oxford, UK 



4 May 2009 



> 

m 

O ■ 1 INTRODUCTION 



in 
o 
o\ 
o 



X 



ABSTRACT 

We propose a model of starburst-driven galactic outflows whose dynamics depends 
on both radiation and thermal pressure. Standard models of thermal pressure-driven 
winds fail to explain some key observations of outflows at low and high redshift galax- 
ies. We discuss a scenario in which radiation pressure from massive stars in a stellar 
population drive a shell of gas and dust. Subsequent supernova (SN) explosions from 
the most massive stars in this population then drive stellar ejecta outward in a rar- 
efied medium, making it collide with the radiation pressure driven shell. The collision 
imparts renewed momentum to the shell, and the resulting re-acceleration makes the 
shell Rayleigh- Taylor unstable, fragmenting the shell. We show that the speed of these 
ballistic fragments can explain some recently observed correlations in Lyman break 
galaxies between wind speed, reddening and star formation rate. 
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Removal of gas from galaxies in the form of outflows con- 
stitute an important phase in the evolution of galaxies, and 
an important mechanism of feedback from galaxies to the 
surrounding medium. At low redshift, outflows have been 
directly observed in many galaxies, in optical to X-rays, and 
at high redshifts,spectroscopic studies have found many ex- 
amples of outflows. Outflows have been ascribed to either 
the stellar population in a galaxy or active galactic nuclei, 
or both (see Veilleux, Cecil, & Bland-Hawthorn 2005, for a 
recent review). 

In the standard scenario of starburst-driven galac- 
tic outflows, thermal pressure builds up in the interstellar 
medium (ISM) of a galaxy owing to supernovae shocks that 
results from formation of stars, and the pressure of the hot 
gas then removes a fraction of the ISM in the form of an 
outflow (Veilleux, Gerald, Bland-Hawthorn 2005, and refer- 
ences therein) . The pressure of the hot gas is directly related 
to the rate of supernovae, or, equivalently to the star forma- 
tion rate. This scenario, however, has difficulty in explaining 
some key observations at both low and high redshifts. 

Starburst-driven winds observed in nearby galaxies 
through Nal absorption lines show a correlation between the 
maximum speed of clouds embedded in the outflow, and the 
star formation rate (Martin 2005). This can be explained by 
ram pressure-driven clouds embedded in the outflow, except 
for very large outflow speeds (Murray, Quataert & Thomp- 
son 2005) that may require radiation pressure (Veilleux, Ger- 
ald, Bland-Hawthorn 2005). Martin (2005) also showed a 



correlation between maximum cloud speed and the rotation 
speed of host galaxies, which she suggested is an indication 
of dominant contribution from radiation pressure. 

Murray, Quataert & Thompson (2005) have discussed 
the basic dynamics of a radiatively driven shell of gas and 
dust. Using these models, Martin (2005) showed that her 
data was indicative of a feedback mechanism on the star 
formation process, in which the star formation rate was reg- 
ulated by gas removal through radiation pressure. 

At high redshifts, Shapley et al. (2003) studied a set of 
~ 800 Lyman break galaxies (LBG) at z ~ 3, and found 
that the doppler shift speed Av = v em — v a b of the Lyman-a 
line is correlated with the star formation rate (as implied by 
the 1Z luminosity), as well as with the reddening due to dust. 
The outflow speed (which is related to Av) was also found 
to be correlated with the equivalent width of low ionization 
absorption lines. Taken together with the correlation with 
the reddening, this suggested that most of the reddening and 
low ionization absorption occurred in the same region, near 
the clouds responsible for the doppler shift. 

Furthermore, Ferrara and Ricotti (2006) suggested af- 
ter studying the outflow speed (V w ~ 0.5Av) and the SFR 
(M) that it was difficult to explain the relation with ther- 
mal pressure-driven winds. The observed relation is approx- 
imately V w oc M 2 / 3 . On the contrary, in the standard SN- 
driven wind via the Sedov- Taylor phase of SN shells, one 
expects V w tx E 1 / 5 , where E oc M is the rate of energy 
deposited by the SN. 

The correlation between reddening and wind speed ob- 
served in LBGs motivates us to seek a model where radiation 
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pressure is a dominant feature. In the case of a purely ra- 
diatively driven shell of gas, the speed is proportional to the 
star formation rate (Murray et al. 2005), so it is not clear 
whether one can explain the above mentioned correlations 
with radiation pressure alone. At the same time, it is reason- 
able to assume that if the pressure of continuum radiation 
comes mostly from massive stars (the continuum dust opac- 
ity being larger in the blue than in the red) , then these stars 
would also explode earlier than the rest of the stellar popu- 
lation. It is conceivable that a radiatively driven wind shell 
is pushed by thermal pressure from SN explosions at a later 
epoch. 

Here we study a hybrid model in which first a radiation 
pressure-driven wind is formed, and a shell moves outward. 
Then, after a few million years, supernova explosions of the 
most massive stars will drive their stellar ejecta through a 
comparatively rarefied medium, and will collide with the ra- 
diatively driven shell. The shell will become vulnerable to 
fragmentation at this stage from Rayleigh- Taylor instabil- 
ity, and the resulting clouds will move ballistically through 
the ISM, with the momentum acquired from the collision. 
We show that this model explains the correlation between 
outflow speed and SFR. In addition, in this model, the dust 
in the cloud, responsible for the dynamics of the radiatively 
driven shell, can explain the relation between reddening and 
outflow speed, as observed in LBGs. 



2 WIND DYNAMICS 

2.1 Phase 1: radiation pressure— driven wind 

We first study a radiatively driven shell of gas and dust. 
We assume a point source of radiation, of luminosity L, and 
assume a spherically symmetric shell of gas at distance r. 
Suppose the total mass (including dark matter) contained 
within this distance is M(r). Setting k to be the opacity of 
matter (dust plus gas) per unit mass of gas (in cm 2 gm -1 ), 
and ignoring gas pressure inside the cavity of the shell, the 
equation of motion of the gas shell is given by 

dV _ GM{r) | kL W 2 
dt r 2 47rr 2 c r 

In a time-independent case, this equation can be re- 
written for a wind, and solved for the wind speed at a given 
radius r, as (ignoring the last term, which refers to decelera- 
tion due to swept-up mass; see Murray, Quataert & Thomp- 
son 2005), 

Here LEd{r) = 47rG -^ f ( r ) c j s the Eddington luminosity rel- 
evant for the given radius r given the mass distribution 
M(r) and opacity k. The star formation rate correspond- 
ing to L Ed (r) is given by M* ~ L Bd (r)/(0.007c 2 ). For 
a singular isothermal mass distribution with velocity dis- 
persion a (M(r) = 2a 2 r/G), this star formation rate is 
~ 50o- 2 00 (r/lOkpc) M Q yr" 1 , for k ~ 500 cm 2 gm" 1 and 
where a = 200<J2oo km s . 



2.2 Dust opacity 

Many observations of galactic outflows indicate the pres- 
ence of dust grains in the outflowing gas. Martin (2005) has 
noted that the absence of Nal emission indicates absorption 
by dust, and that cold dust is often found in starburst galax- 
ies at larger scale heights than those of the stars. Opacities 
of order k, ~ 500 cm 2 gm -1 can be provided by dust grains 
through continuum absorption and scattering (Draine & Lee 
1984). The absorption opacity of gas mixed with dust is ex- 
pected to be larger at shorter wavelengths: k ~ 600-100 
cm 2 gm" 1 at A ~ 1000-6000 A (Li & Draine 2001). Obser- 
vations of reddening of local (Calzetti 2001 and references 
therein) and high redshift starburst galaxies (Sanders & 
Mirabel 1996), including LBGs (Adelberger & Steidel 2000) 
also suggest dust opacities of this order. In our calculations 
described below, we assume dust opacity of order k ~ 500 

2 —1 

cm g . 

2.3 Evolution in luminosity 

Consider the case of a single burst of star formation with 
an initial mass function (IMF) <j>(m) = 4>\(jn/m\)~ < ^ +x \ 
which is normalized as M a = M J m<f>(m)dm, where M is 
the total mass of the stellar population. The total luminosity 
of main sequence stars in a stellar population evolves as, 

L(t) oc / <t>(m)l(m)dm ' w (t/n) t 1 -") , (3) 

/ a — x 

where Z(m) = li(m/m\) a is the main sequence luminosity 
of a star with mass m, and a ~ 3.5. The main sequence 
lifetime of a star with mass mi is denoted here by t\. Stars 
on the giant branch contribute more to the luminosity, and 
their contribution can be incorporated by including a mul- 
tiplicative factor [1 + G(t)], where G{t) = ggf ; ( '^ )t , and 
l g , Tg are the giant branch luminosity and lifetime on the gi- 
ant branch, and mtn is the turn-off mass at time t (Tinsley 
1980). The ratio G(t) is of order unity for short wavelengths 
(although larger for longer wavelengths, where giant branch 
stars contribute most), and so we multiply the expression in 
equation |3] by a factor 2. 

In the case of continuous star formation, with a constant 
SFR, it can be shown that the luminosity evolves as 

w [l + G(t)]M4> imi h (Vn) l5f ^ 4 xl0 9 L s M 1 tg- 14 ,(4) 

a — x 

where t a = t/l Myr, Mi = M/l M yr -1 , and where we 
have assumed a Salpeter IMF (x = 1.35), and that G(t) ~ 1. 
We use this expression for luminosity in equation [1] for our 
calculations below. 

2.4 Phase 2: supernova explosions 

Next, we consider the evolution of the stellar ejecta from 
supernova explosions of the most massive stars in this pop- 
ulation. The ejecta will move through a rarefied medium in 
the interior of the radiatively driven shell, since a large frac- 
tion of gas is expected to be swept out earlier by radiation 
pressure. We consider supernovae resulting from the early 
stages of a starburst, within an interval t sn ~ 3 Myr, which 
is the main sequence lifetime of stars with M > 15M©. 
Assuming a Salpeter IMF, one gets a total mass of ejecta 
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(ignoring the remnant masses), as M e j ~ 0.1Mt 3n , consid- 
ering supernovae for M > 15M©, and a lower and upper 
mass cut-offs of O.f and 125 Mq. The total energy output is 
related to the number of SN, ~ 0.003 Mt sn , again assuming 
a Salpeter IMF and supernovae for M > 15Mq . We assume 
an energy output of 10 51 ergs per supernova, so that the 
total energy output is E e j = 10 51 erg x 0.003 Mt sn . 

We also assume that the evacuated and rarefied ISM in 
the interior of the radiatively driven shell has a particle den- 
sity n = n r . In the Sedov- Taylor phase (which is valid since 
the ejecta travel with high speed until colliding with the ra- 
diatively driven shell, and since the shell does not sweep up 
much matter in the rarefied medium to slow it down), we 
have for the evolution of the outer shell of the ejecta (e.g., 
Shull & Silk 1979), 

R ej « 9.89 x 10" 2 kpc E° e 2 bl n"?' 2 tg 4 
V ej « 38.69 kms -1 £°jf 51 n^?' 2 tg °' 6 , (5) 
where Eej, 51 = E e j /10 51 erg, and n rt i = n r j\ cm -3 . 



2.5 Re-acceleration of the shell and growth of 
instability 

The radiatively driven shell is first accelerated for a short 
period, and then soon begins to decelerate as the shell radius 
r increases and the radiation pressure (which scales as r~ 2 ) 
falls. At the same time, the shell of ejecta from subsequent 
SN events moves swiftly through the rarefied medium inside 
the radiatively driven shell, and catches up with it at a time 
tcoii, when they collide. The speed of the ejecta is in general 
much larger than the shell speed, and although the ejecta 
mass is very small, the momentum imparted to the shell by 
the colliding ejecta is large. Therefore, the shell is carried 
forward by the increased momentum from the ejecta coming 
up behind it. The resulting speed of the shell is given by 

V = {MejVsn + M shell V sheH )/(M ej + M sheH ) , (6) 

where M s hell and V s hell are the swept-up mass in the shell 
and the speed of the shell at the time of collision, and V sn 
is the speed of the ejecta at that time. 

After the collision, the wind shell is re-accelerated, and 
we have the case of a heavy fluid (in the shell) above a rar- 
efied medium, both being acted upon by an effective gravity 
resulting from the acceleration of the shell. The shell ma- 
terial is therefore subject to Rayleigh- Taylor (RT) instabil- 
ities, and is liable to fragmentation. 

Even before the collision, the expanding shell is prone 
to dynamical and gravitational instabilities. Vishniac (1983) 
estimated the shortest growth time for gravitational insta- 
bility to be given by t g ~ c s /irGa, where a is the surface 
mass density. A lower limit on the growth time for instability 
is obtained by assuming a cool shell with T ~ 10 4 K. 

After the collision, the shell is re-accelerated, and be- 
comes Rayleigh- Taylor unstable. The effective gravitational 
acceleration at this time is g B jj ~ v 2 /r. The time scale of 
Rayleigh- Taylor instability for clouds of size A = 2n/k is 
given by, 

trt ~ {g B ffk)~ 1/2 ~ (X/2-Kg eff ) 1/2 , (7) 
where k is the wavenumber. 



3 RESULTS 

We follow the evolution of the radiatively driven shell (eqn 
[TJ numerically and track the trajectory of SN ejecta (eqn[5]), 
and determine the time of collision (t co u) between the two. 
For the dark matter profile, we have assumed the galaxies 
to be of total mass M = 10 12 Mq, at z = 3, and assumed 
a Navarro-Frenk- White halo profile, with a concentration 
parameter c = 3 (Navarro, Frenk, White 1996). Our as- 
sumption is motivated by the estimate made by Adelberger 
et al. (2005) from the clustering of LBGs and comparison 
with N-body simulations. The median halo mass of LBGs 
at z ~ 3 is estimated to be io u - 86±0 - 3 Mq. Although there 
is considerable uncertainty in the dynamical mass estimate 
of LBGs (see Weatherley & Warren 2005), our results are 
not sensitive to the gravitational deceleration of the shell by 
dark matter, and therefore, to the assumption of total mass 
and its profile. 

Figure 1 shows the evolution of the shell radius (solid 
line, top panel) and speed (solid line, bottom panel) with 
time, for a star formation rate of M = 100 Mq yr _1 , a 
continuum dust opacity n = 500 cm 2 gm . The particle 
density in the undisturbed ISM (that the radiatively driven 
shell sweeps through) is assumed to be n g = 5 cm -3 , and 
the density inside the shell is assumed to be n r = 0.01 cm' 3 . 
It is seen that the shell accelerates for an initial brief period 
(<C 1 Myr), after which it continuously decelerates, mostly 
owing to the decrease in radiation pressure with distance, 
and to some extent, to the growing weight of the swept-up 
mass in the shell. 

The dotted lines in Figure 1 show the corresponding 
evolution for the stellar ejecta, which we assumed to expand 
outward after a period of t sn ~ 3 Myr. The rarefied medium 
inside the radiatively driven shell helps the ejecta move with 
high speed, and catch up with the shell within a short time. 

Figure 2 plots the (logarithm of the) ratio of the short- 
est time-scale for instability to the age of the shell, as a 
function of time. The dotted line shows the shortest time- 
scale for gravitational instability, assuming a cool shell with 
gas at T ~ 10 4 K. The surface mass density is calculated 
from the swept-up mass within r, the shell radius. The fig- 
ure shows that gravitationally unstable modes takes ^> 10 
Myr to grow, although the time-scale for growth steadily 
decreases with time, due to the increasing surface mass den- 
sity. 

After the collision, the shortest time-scale for RT in- 
stability is calculated by using equation [7] for instabilities of 
wavelength ~ 2nr, where r is the shell radius. This is the 
largest wavelength for instabilities in the shell, and therefore 
gives the time-scale for the fastest growing mode. This time- 
scale at the time of collision is shown as a point. The effective 
gravitational acceleration at this time is g e ff ~ v 2 jr ~ 10 -7 
cm s -2 , for M ~ 100 Mq yr _1 , making the shortest time- 
scale for instability growth comparable to the age of the 
shell (a few Myr). The re-acceleration of the shell, there- 
fore, makes it prone to RT instabilities. This is likely to feed 
gravitational instabilities on scales that are gravitationally 
stable at this stage. 

The size of the resulting fragments will be 1-1.5 kpc, 
the shell radius at this point of time, as plotted in Figure 
3 as a function of the outflow speed for the assumed val- 
ues of the parameters. We then identify the speed of the 
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Figure 1. The top panel plots the shell radius (kpc) with time 
(Myr) (solid line), for a star formation rate of 100 Mq per year, 
with ambient gas density n g = 5 per cc, and for opacity k = 500 
cm 2 gm _1 . The dotted line shows the evolution of the hot gas 
bubble, assuming a density inside the shell n r = 0.01 per cc. The 
bottom panel shows the corresponding speed of the shell (solid 
line) and the hot bubble (dotted line). 




Figure 2. The ratio between the instability growth time-scale 
and age of the shell (in logarithmic units) is plotted with time. 
The dotted line shows the shortest time-scale for gravitational 
instability, assuming a shell temperature of 10 4 K. For most of 
the time span, the shell decelerates and so there is no RT insta- 
bility, but after the collision, it is accelerated, and the ratio of RT 
timescale and age of the shell is shown as a point at the time of 
collision. 
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Figure 3. The shell radius at t co u as a function of outflow speed 
(in km/s) for the same set of parameters as in Figure 1. 

cloud after collision as the outflow speed, V w . We can then 
compare our results with data from Shapley et al. (2003), 
for different assumptions of parameter values. The results 
of the present model are insensitive to the launching radius 
and initial speed of the shell, and depends on the assumed 
values of dust opacity (ft), the densities of gas outside (n g ) 
and inside (n r ) the shell. We found that the best case sce- 



nario needs k 



500 cm 2 gm" 1 , 



5 cm ' , and n r 



0.01 



The top panel of Figure 4 plots the resulting outflow 
speed with SFR for these parameter values , where we show 
the result of our calculation with a solid curve, along with 
data points from Shapley et al. (2003) (their Table 3, with 
h — 0.7). The dashed line in the top panel shows the slope 
of the expected relation between outflow speed and SFR 
(scaled arbitrarily) in the case of a shell driven by hot gas 
(V w ocS 1 / 5 ). 

Another natural prediction of this scenario is a correla- 
tion between outflow speed and reddening due to dust ab- 
sorption. The outflow speed is determined by the the com- 
bined momentum of the radiative shell and the ejecta, which 
is related to the total mass in the shell, and therefore, to 
the amount of dust in the shell. We can quantify this cor- 
relation by estimating the amount of reddening using the 
column density of material in the shell. The bottom panel 
of Figure 4 shows the expected reddening (E(B-V)) as a 
function of the cloud speed for three different choices of the 
reddening versus column density relation. The solid line at 
the bottom of the figure uses the relation between hydrogen 
column density and reddening as found in the Small Magel- 
lanic Cloud (SMC) (N H /E(B - V") = 4.4 ± 0.7 x 10 22 cm 2 
mag -1 ; Bouchet et al 1985), with the dotted lines showing 
the 1 — a spread. The short-dashed line uses the calibra- 
tion from observations of the Large Magellanic Cloud (LMC; 
N H /E{B-V) = 2±0.5 x 10 22 cm 2 mag" 1 ; Koornneef et al 
1982), and the long-dashed line uses the Galactic ISM data 
(N H /E(B - V) = 4.93 ± 0.28 x 10 21 cm 2 mag" 1 ; Diplas & 
Savage 1994). 



4 DISCUSSION 

The top panel of Figure 4 shows that the scenario described 
here, along with suitable assumptions of the values of key 
parameters (n = 500 cm 2 gm -1 , n g — 5 cm~ 3 , and n r — 0.01 
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Figure 4. The top panel plots the outflow speed with SFR with 
a solid curve. The dashed curve is the expected result from a 
thermal pressure-driven wind. The bottom panel plots the red- 
dening with outflow speed using calibration of column density 
and E(B-V) from SMC (solid line), LMC (short-dashed line) and 
Galactic ISM (long-dashed line). The dotted lines show the 1-cr 
error margin for the calibration for SMC type dust. 



cm ), can explain the observations of outflows speed and 
SFR observed in LBGs. 

While the calibration of column density and reddening 
from SMC data explains the data of reddening and out- 
flow speed in LBGs from Shapley et al. (2003), the expected 
reddening from Galactic ISM calibration is much larger. It 
should be kept in mind that dust grains are destroyed in 
ISM shocks by a variety of processes, from energetic grain- 
grain collisions and from collisions between gas ions and dust 
grains. Calculations show that large grains (with sizes larger 
than ~ 0.05 fx) are destroyed in shocks with speed larger 
than ~ 200 km/s (Jones et al. 1997). In the present case, 
dust grains embedded in the radiatively driven shell coast 
with a speed ~ 200 km s _1 for a few Myr before being hit by 
stellar ejecta (see bottom panel of Figure 1). Afterwards, the 
fragmented clouds traverse the rest of the ISM at high speed 
(v ~ 300-400 km s _1 ) and leave the host galaxy within a 
few Myr, because the typical half-light radius of LBGs is of 
order ~ 1 kpc (see below), comparable to the distance where 
clouds form in the present scenario. Dust grains in the cloud 
will therefore be destroyed to a certain extent, especially the 
larger grains, and the dashed and dotted lines in Figure 4 
(bottom panel), corresponding to Galactic and LMC-type 
dust, will come closer to the LBG data. 

Another interesting possibility is that due to the de- 
struction of large grains, the grain size distribution in these 
clouds will be dominated by smaller (< 500 A) grains, and so 



the total surface area of grains will increase, thereby increas- 
ing the reddening for shorter wavelength ultraviolet contin- 
uum. The SMC-type extinction curve shows more redden- 
ing at wavelengths shorter than the 2175A feature than the 
Milky Way extinction curve (e.g., Pei 1992), and so the solid 
line in Figure 4, corresponding to SMC-type dust, is likely to 
show increased reddening for larger outflow speeds, thereby 
bringing the predicted curve closer to the data points in 
Shapley et al. (2003). At the same time, the smaller grain 
size would show less reddening in the optical. Therefore, the 
frequency dependence in reddening for outflows in LBGs will 
likely change with redshift, as the rest frame wavelength for 
observed optical radiation changes with redshift. 

An important caveat arises from the fact that in com- 
paring with the data, we have identified the outflow speed 
with AV/2, where AV = \V em - V abs \. Adelberger et al. 
(2005) have shown that there are asymmetries between 
Vem - V ne b and V a bs - V neb , where V neb is the speed of the 
ISM, indicating that outflow speed can significantly differ 
from AV/2, which we have assumed for simplicity. 

The scenario discussed here has several implications. 
Firstly, the shell fragments are generated at radii ~ 11.5 
kpc (Figure 3), comparable to the half-light radii of LBGs. 
Bouwens et al. (2004) estimated the mean half-light size 
of LBGs at z ~ 3 as ~ 1.3 kpc. Therefore, the fragments 
are likely to penetrate the rest of the galactic ISM and de- 
posit metals into the intergalactic medium (IGM), instead 
of falling back on to the host galaxies. Galactic outflows 
have been invoked to explain the metallicity of low column 
density Lyman-a absorbers in the IGM, and several authors 
have expressed concern that this process of metal enrichment 
through outflows might (1) increase the number of absorp- 
tion lines more than necessary for consistency with models 
of structure formation, and (2) dynamically disturb the ab- 
sorption systems. It has been suggested that these problems 
could be overcome if the wind shell fragments, so that the 
fragments can coast to large distances (e.g., Theuns, Mo & 
Schaye 2001). (Recent detailed numerical studies, however, 
find that galactic winds may not pose a problem for the 
interpretation of Lyman-a absorption systems as was pre- 
viously thought (Cen et al. 2005).) We note that galactic 
wind-driven shells in the scenario discussed here undergo 
fragmentation even while traversing the ISM, and the frag- 
mented shell is likely to spread metals in the IGM in an 
inhomogeneous manner, and not dynamically disturb the 
gaseous filaments responsible for Ly-a absorption systems. 
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